Thermodynamics of condensed matter with strong pressure-energy correlations 
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We show that for any hquid or sohd with strong correlation between its NVT virial and potential- 
energy equilibrium fluctuations, the temperature is a product of a function of excess entropy per 
particle and a function of density, T = f{s)h{p). This implies that 1) the system's isomorphs (curves 
in the phase diagram of invariant structure and dynamics) are described by h{p)/T = Const., 2) 
the density-scaling exponent is a function of density only, 3) a Griineisen-type equation of state 
applies for the configurational degrees of freedom. For strongly correlating atomic systems one 
has h{p) = ^^ CnP^ '^ in which the only non-zero terms are those appearing in the pair potential 



expanded as v(r) = 5]]„ '^n 
confirm the theory. 



Molecular dynamics simulations of Lennard- Jones type systems 



The class of strongly correlating liquids was introduced 
in Refs. [l| and y- These liquids are defined by hav- 
ing a correlation coefficient above 0.9 of the constant- 
volume equilibrium fluctuations of virial W and poten- 
tial energy U . The WU correlation coefficient varies with 
state point, but we found from computer simulations 
that a system has either poor WU correlations in the 
entire phase diagram or is strongly correlating at most 
of its condensed-phase state points X"^ ■ ^^'^ der Waals 
and metallic liquids are generally strongly correlating, 
whereas hydrogen-bonded, ionic, and covalently bonded 
liquids are generally not. The solid phase is at least as 
correlating as the liquid phase. Theoretical arguments, 
numerical evidence, and experiments show that strongly 
correlating liquids are simpler than liquids in general [l|- 

The simplicity of strongly correlating liquids compared 
to liquids in general derives from the fact that the for- 
mer have "isomorphs" in their phase diagram. Two state 
points with density and temperature (pi, T\) and (p2, T2) 
are termed isomorphic [3| if all pairs of physically relevant 
microconfigurations of the state points that trivially scale 
into one another (i.e., p^ r,^ — pj ^l for all particles 
i), have proportional configurational Boltzmann factors: 
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Only inverse-power law liquids ^] have exact isomorphs 
(here C12 — 1), but as shown in Appendix A of Ref. 
y a system is strongly correlating if and only if it has 
isomorphs to a good approximation. 

The invariance of the canonical probabilities of scaled 
microconfigurations along an isomorph has several im- 
plications, for instance [i|-y|' 1) The excess entropy and 
the isochoric specific heat are isomorph invariants, 2) 
the reduced-unit dynamics is isomorph invariant for both 
Newtonian and stochastic dynamics, 3) all reduced-unit 



static correlation functions are isomorph invariant, 4) 
a jump between isomorphic state points takes the sys- 
tem instantaneously to equilibrium. Using reduced units 
means measuring length in terms the unit p^^'^ where 
p = N/V is the particle density, and time in units of 
p^^'^y^m/kBT where m is the average particle mass. 
Since isomorphs are generally approximate, isomorph 
properties are likewise rarely rigorously obeyed. 

All thermodynamic quantities considered below are ex- 
cess quantities, i.e., in excess of those of an ideal gas at 
the same density and temperature. Thus S is the excess 
entropy (5 < 0), Cy is the excess isochoric specific heat, 
p is the excess pressure (i.e., p — W/V), etc. 

Briefly, the reason that S and C'v are isomorph invari- 
ants is the following [3[ . The entropy is determined by the 
canonical probabilities, which are identical for scaled mi- 
croconfigurations of two isomorphic state points. From 
Einstein's formula Cy = {{AU)'^)/kBT^ the isomorph 
invariance of Cv follows easily by taking the logarithm 
of Eq. ([1]) and making use of the isomorph invariance of 
scaled microconfiguration probabilities. 

Since S and Cy are invariant along the same curves 
in the phase diagram, C'v is a function of S: Cv — 
(t)iS). Thus T{dS/dT)v = (t){S) or at constant volume: 
dS/(j){S) = dT/T. Integrating this leads to an expres- 
sion of the form ip{S) = ln(T) -|- k(V), which implies 
T = exp[f/'(S')] exp[— A:(T^)]. The generic version of this 
involves only intensive quantities (s = S/N): 



T = f{s)h{p) . 



(2) 
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For inverse power law interactions (ex r^") the entropy 
is well-known to be a function of p^ /T where 7 — n/3: 
S = K^p^ /T). Applying the inverse of the function K , 
shows that these perfectly correlating systems obey Eq. 
dll) with h{p) = p^. 

The thermodynamic separation identity Eq. ^ is the 
main result of this paper. We proceed to discuss some 
consequences and numerical tests. 

1. Density scaling 
Since entropy is an isomorph invariant, it follows from 
Eq. ^ that the variable characterizing an isomorph 



may be chosen as h{p)/T. In particular, the reduced 
relaxation time f, which is also isomorph invariant, may 
be written for some function G 
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This is the form of "density scaling" proposed by Alba- 
Simionesco et al. in 2004 from different arguments [10|; 
at the same time Dreyfus et al., as well as Casalini and 
Roland, favored the more specific form f = G{p'^ /T) \\^ . 
Isochrones for many supercooled liquids and polymers 
follow to a good approximation the latter "power-law 
density scaling" relation [il[. For large density changes, 
however, it was recently shown that the density-scaling 
exponent generally varies in both simulations and exper- 
iment [131; these cases conform to the more general Eq. 
©. 

2. An expression for the density- scaling exponent 
The general, state-point dependent density-scaling expo- 
nent 7 is defined 0, 3 by 
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The physical interpretation of Eq. ^ is the following. 
If density is increased by 1%, temperature should be 
increased by 7% for the system to have the same en- 
tropy and reduced relaxation time. Equation ([2]) implies 
dlnT = dln/(s) + dln/i(p); thus along an isomorph - 
where s and f are both constant ~ one has d In T = d In ft,. 
Via Eq. @ this implies 



U == F{S)h(p) + k{p) leads to Eq. ©, in which 7(p) is 
given by Eq. ([5]). 

^. The isoniorphs of atomic systems 
We consider now predictions for systems of "atomic" 
particles interacting via pair potentials of the form [1^ 
(where r is the distance between two particles) 
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For simplicity of notation the case of identical particles 
is considered, but the arguments generalize trivially to 
multicomponent systems. Consider the thermal average 
(r^"). Switching to reduced units defined by f = p^^^r, 
we have (r~") = (f~")p"/'^. Since structure is isomorph 
invariant in reduced units, (f^"') is an isomorph invari- 
ant. Consequently, it is a function of any other isomorph 
invariant, for instance the entropy: {r~"') = GniS). 
Noting that the average potential energy is a sum of 
Eq. ([7]) over all particle pairs, we conclude that (where 
i/„(^)cx«„G„(5)) 
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Taking the derivative of this equation with respect to 
temperature at constant volume leads to 
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(5) The left hand side is T {dS/dT)y, so Eq. © implies 



In particular, 7 depends only on density: 7 = ^{p) [3]. 

3. Configurational Griineisen equation of state 
The Griineisen equation of state expresses that pressure 
equals a density-dependent number times energy plus a 
term that is a function of density only J^] . This equa- 
tion of state is used routinely for describing, in particu- 
lar, solids under high pressure. We proceed to show that 
strongly correlating systems obey the configurational ver- 
sion of the Griineisen equation of state, which as sug- 
gested by Casalini et al. [1J| has the density-scaling ex- 
ponent as the proportionality constant [3, [j] 



W = -i{p)U + $(p) , 
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To prove this, note first that {dU/dS)^ = T = f{S)h{p) 
by integration implies U — F{S)h{p) + k{p) where 
F'{S) = f{S) {S is the extensive entropy). Since W = 
{dU / din p) g (which follows from the standard identity 
TdS = dU+pdV), we get W = F(S')dft/dlnp-|-dfc/dlnp. 
Substituting into the latter expression F{S) isolated from 



T = ^if;(5)/ 
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This is consistent with the thermodynamic separation 
identity Eq. ([2|) only if all the functions H'^{S) are 
proportional to some function, i.e., if one can write 
H[J^S) = C„(/)(S'). We identify (/)(S') as the function /(s) 
of Eq. ([2]), which means that 



h{p) = ^C„p"/3 



(11) 



Thus for strongly correlating atomic liquids, the thermo- 
dynamic function h{p) has an analytical structure, which 
is inherited from v{r) in the sense that the only non-zero 
terms of h{p) are those corresponding to the non-zero 
terms of v{r). Note that not all systems with potentials 
of the form Eq. ([7]) are strongly correlating and that the 
derivation applies only if this is the case. 
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99.9% in its entire phase diagram. At low densities 
{p <C 1) the repulsive LJ fluid behaves as an r~^ fluid, 
whereas it for p ::^ 1 is effectively an r^^^ fluid. Thus 
the density-scaling exponent "f{p) varies from 2 to 4 as 
density increases, a much larger variation than that of 
previously studied strongly correlating systems. 



Since h(p) is only defined within an overall multiplica- 
tive constant, one can write for the repulsive LJ fluid 
h{p) = ap'^ + {1 - a)p'^. This leads via Eq. ^ to 
7o = 2 + 2a, implying that 



FIG. 1: The thermal average of r~^^ versus that of — r~® in 
reduced units for a large range of state points of the Kob- 
Andersen binary Lennard-Jones liquid simulated with 1000 
particles (saa = (JAA = 1). These quantities correspond 
to Hi2{S) and HsiS) in Eq. ((Hjl. The theory predicts that 
H'i2{S) oc H(i{S), implying that all data points should fall 
onto a common line according to Hi2{S) = aHaiS) + /3. 

As an illustration we present results from simulations 
of the Kob-Andersen binary Lennard-Jones (KABLJ) liq- 
uid [l6f|, which is strongly correlating at its condensed- 
phase state points [i|-|j] . The application of the above to 
LJ systems predicts that H[2{S) oc iJg(S'), where Hi2{S) 
is the reduced coordinate average of the r~^^ term of U, 
etc. Lrtegrating this leads to Hi2{S) — aHQ{S) + (3, im- 
plying that if the repulsive term in U is plotted against 
the attractive term in reduced units, all points should 
fall onto a common line. Figure [T] presents data where 
density was changed by a factor of eight and tempera- 
ture a factor of 40,000. The data collapse is good but 
not exact, which reminds us that the relations derived 
are approximate. 

The theory implies a simple mathematical description 
of the isomorphs in the {p, T) phase diagram. From the 
fact that the potential energy contains only r~^'^ and r~^ 
terms, it follows that h{p) — Ap'^ — Bp^. Consequently, 
LJ isomorphs are given by 



Ap'^ - Bp^ 
T 



= Const. 



(12) 



The invariance of the Boltzmann statistical weights of 
scaled microconfigurations implies that an isomorph can- 
not cross the liquid-solid coexistence curve. In particular, 
the coexistence curve is itself predicted to be an isomorph 
[3| , which was recently confirmed by simulations of gen- 
eralized L J liquids [J, [T^l . Consequently the coexistence 



line is given by Eq. (|T2)) . This validates a recent conjec- 
ture of Khrapak and Morfill [ll]. 

5. Predictions for the repulsive Lennear- Jones fluid 
As a final illustration we consider the "repulsive" single- 
component LJ fluid defined by the pair potential {r~^'^ + 
r~®)/2, a system with WU correlation coefficient above 



h{p) = (7o/2 - 1)/ + (2 - 7o/2)p2 



(13) 



Our simulations identified from the expression 70 — 
(AVFAC/)/((At/)2) g the exponent 70 = 3.56 at the 
state point (p, T) = (1, 1). Equation p^ with 70 = 3.56 
was tested in two different ways. First, we compared at 
each state point along an isomorph the exponent 7(p) 
predicted from Eqs. ([5]) and P^ with that calculated 
from the fluctuations via 7 = (AW AU) / {{AU)'^) (right 
panel of Fig. 2). The left panel presents a second test 
of Eq. (|13p by showing results from simulating five tem- 
peratures at p = 1, plotting for each temperature in- 
stantaneous values of the potential energy versus the po- 
tential energy of the same microconfigurations scaled to 
three other densities (p = 0.5,1.6,2.0). The theory be- 
hind the observed straight lines is the following. Consider 
two isomorphic state points (poj^o) and {p,T) and sup- 
pose each temperature is changed a little, keeping both 
densities constant. If the two new state points are also 
isomorphic, the entropy change is the same for both: 
dUo/To = dU/T. This implies dU/dUa = T/Tq, i.e., 
{dU/dUQ)p^p = T/To. Since h{p)/T is constant along 
an isomorph, this implies {dU/dUo) — h{p)/h{po)- 
Integrating this at constant po and p leads to U = 
[h{p) / h{pQ)]Uo + 4>{po,p). In our case po = 1 and 
h^po) = 1. Thus plotting U versus Uq is predicted to 
result in straight lines with slope h{p) (yellow asterices 
in the left panel of Fig. ^. The scaled state points are 
isomorphic to the original p = 1 state points, with tem- 
peratures given by T = Toh{p). Via the "direct isomorph 
check" [3| this implies that the scaled microconfigurations 
form elongated ovals with slope h{p). 
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In summary, we have shown that for strongly correlat- 
ing liquids or solids, temperature separates into a func- 
tion of entropy times a function of density. For these 
systems the energy scale is consequently determined by 
density alone. It is an open question whether, conversely, 
the thermodynamic separation identity Eq. ([2]) implies 
that the system in question is strongly correlating. 



FIG. 2: "Multiple direct isomorph check" applied to simula- 
tions of TV = 1000 particles of the repulsive LJ fluid defined 
by the pair potential (r~^^ -|- r~^)/2. The left panel shows 
the potential energies of pairs of microconfigurations, where 
the potential energy of a given microconfiguration at density 
1.0 is denoted C/(1.00) and that of the same microconfigura- 
tion scaled to density p is denoted U{p) (p = 0.5; 1.6; 2.0). 
This was done for T = 0.6; 0.8; 1.0; 1.2; 1.4. The black 
lines are the predictions (see the text) with slopes deter- 
mined via Eq. 1)13^ from the fluctuations calculated at the 
state point {p,T) = (1,1) marked by an arrow. The right 
panel shows the density-scaling exponent for each state point 
along an isomorph predicted from Eqs. ([5]) and (|13|1 (full 
curve) and the exponent calculated via the fluctuation for- 
mula 7 = {AWAU)/{{AUf) 3 (red crosses). The arrow 
marks the state point {p,T) = (1, 1). 



The centre for viscous liquid dynamics "Glass and 
Time" is sponsored by the Danish National Research 
Foundation (DNRF). 



[1] N. P. Bailey et al, J. Chem. Phys. 129, 184507 (2008); 
N. P. Bailey et al., J. Chem. Phys. 129, 184508 (2008); 
T. B. Schr0der et al, J. Chem. Phys. 131, 234503 (2009). 

[2] U. R. Pedersen et al, Phys. Rev. Lett. 100, 015701 
(2008); N. Gnan et al, Phys. Rev. Lett. 104, 125902 
(2010); U. R. Pedersen et al, Phys. Rev. Lett. 105, 
157801 (2010). 

[3] N. Gnan et al, J. Chem. Phys. 131, 234504 (2009). 

[4] T. B. Schr0der et al, J. Chem. Phys. 134, 164505 (2011). 

[5] U. R. Pedersen et al, J. Non-Cryst. Solids 357, 320 
(2011). 

[6] D. Gundermann et al. Nature Physics 7, 816 (2011). 

[7] T. S. Ingebrigtsen, T. B. Schr0der, and Jeppe C. Dyre, 
'arXrv:1111.3557 (2011). 

[8] O. Hirschfelder, C. F. Curtiss, and R. B. Bird, Molec- 
ular theory of gases and liquids (Wiley, New York, 
1954); J. P. Boon and S. Yip, Molecular hydrodynam- 
ics (McGraw-HiU, New York, 1980); J. S. Rowlinson and 
B. Widom, Molecular theory of capillarity (Clarendon, 
Oxford, 1982); M.P. Allen and D.J. Tildesley, Computer 
simulation of liquids (Oxford Science Publications, Ox- 
ford, 1987); D. Chandler, Introduction to modern sta- 
tistical mechanics (Oxford University Press, New York, 
1987); P. G. Debenedetti, Metastable liquids: Concepts 
and principles (Princeton University Press, Princeton, 



NJ, 1996); N. H. March and M. P. Tosi, Introduction 
to liquid state physics (World Scientiflc Publishing, Sin- 
gapore, 2002); J.-L. Barrat and J. -P. Hansen, Basic con- 
cepts for simple and complex liquids (Cambridge Univer- 
sity Press, Cambridge, England, 2003); J. -P. Hansen and 
J. R. McDonald, Theory of simple liquids, 3rd ed. (Aca- 
demic, New York, 2005). 
[9] O. Klein, Medd. Vetenskapsakad. Nobelinst. 5, 1 (1919); 
T. H. Berfin and E. W. MontroU, J. Chem. Phys. 20, 75 
(1952); W. G. Hoover et al, J. Chem. Phys. 52, 4931 
(1970); W. G. Hoover, S. G. Gray, and K. W. John- 
son, J. Chem. Phys. 55, 1128 (1971); Y. Hiwatari et al. 
Prog. Theor. Phys. 52, 1105 (1974); D. M. Heyes and 
A. C. Branka, J. Chem. Phys. 122, 234504 (2005); A. 
C. Branka and D. M. Heyes, Phys. Rev. E 74, 031202 
(2006). 

[10] C. Alba-Simionesco, D. Kivelson, and G. Tarjus, J. 
Chem. Phys. 116, 5033 (2002); C. Dreyfus et al, Phys. 
Rev. E 68, 011204 (2003); C. Alba-Simionesco et al, 
Europhys. Lett. 68, 58 (2004); R. Casalini and C. M. 
Roland, Phys. Rev. E 69, 062501 (2004). 

[11] C. M. Roland et al. Rep. Prog. Phys. 68, 1405 (2005); 
G. Floudas, M. Paluch, A. Grzybowski, and K. L. Ngai, 
Molecular Dynamics of Class-Forming Systems: Effects 
of Pressure (Advances in Dielectrics, Springer, 2010); 



D. Fragiadakis and C. M. Roland, J. Chem. Phys. 134, 

044504 (2011). 

[12] L. B0hling et tt/., larXiv:1112.1602K 201lV 

[13] M. Born and K. Huang, Dynamical Theory of Crystal 

Lattices (Oxford University Press, Oxford U.K., 1954); 

M. Ross and D. A. Young, Annu. Rev. Phys. Chem. 44, 

61 (1993); L. Buralcovsky and D. L. Preston, J. Phys. 

Chem. Solids 65, 1581 (2004). 
[14] R. Casalini, U. Mohanty, and C. M. Roland, J. Chem. 



Phys. 125, 014505 (2006). 
[15] Y. Rosenfeld, Phys. Rev. A 26, 3633 (1982). 
[16] W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 

(1994). 
[17] A. Ahmed and R. J. Sadus, J. Chem. Phys. 131, 174504 

(2009). 
[18] S. A. Khrapak and G. E. MorfiU, J. Chem. Phys. 134, 

094108 (2011). 



